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Abstract 

We present the results of an analytic next -to-next -to leading order QCD calculation of the 
non-singlet anomalous dimension 7^g (iV) and the coefficient functions C^zX-ZV) associated 
to the deeply inelastic structure functions F2(x, Q 2 ) and Fl(x, Q 2 ) for the Mellin moment 
N = 16. Comparisons are made with results in the literature. 



1 Introduction 



One of the most precise methods to determine the strong coupling constant a s (Q 2 ) is to measure 
it from the scaling violations of the structure functions in deeply inelastic scattering. The present 
level of experimental accuracy is of order 5a s (M|) ~ ±0.001 [1]. Since the theory errors due to 
factorization and renormalization scale uncertainties are of 5a s (M^) sc ~ ±0.005 [2] in next-to- 
leading order the knowledge of the 3-loop anomalous dimensions is required to reduce this error 
to the level reached by experiment. 1 

The one-loop anomalous dimensions were calculated in [5] and the one-loop coefficient func- 
tions were obtained in a series of subsequent calculations, see Ref. [6] for a summary of the 
results. Later the two-loop anomalous dimensions [7] and the two-loop coefficient functions [8] 
were evaluated applying different methods. The calculation of the anomalous dimensions and 
Wilson coefficients in three loop order is a very difficult task. The first calculations were per- 
formed for a series of sum rules and the individual moments for N = 2, 4, 6, 8 [9] in the non-singlet 
and singlet case as well as for N = 10 in the non-singlet case using the symbolic manipulation 
program FORM [10]. With growing computational power the moments N = 10, 12 in the sin- 
glet case and N = 12, 14 in the non-singlet case as well as the moments N = 3, 5, 7, 9, 11, 13 
associated to the structure function xF 3 (x, Q 2 ) could be calculated [11]. Very recently, the 3- 
loop non-singlet [12] and singlet [13] anomalous dimensions were calculated for deep-inelastic 
scattering off unpolarized nucleons. 

In this paper we calculate the 16th moment of the non-singlet splitting function, 
1ns' + (N, Q 2 ), and the Wilson coefficients of the structure functions F 2 (x,Q 2 ) and F L (x,Q 2 ), 
C^(N,Q 2 ), for the case of pure photon exchange to 3-loop order using the MINCER formal- 
ism [14]. This calculation was started before complete 3-loop results became available and 
serves as an independent check of the result for the non-singlet anomalous dimension given in 
Ref. [12] and predicts the 16th moment for the non-singlet coefficient function, [15,16]. The 
paper is organized as follows. In section 2 we give an outline of the basic formalism. The renor- 
malization of the hadronic matrix element is described in section 3. The results of the calculation 
are given in section 4 and section 5 contains the conclusions. 



2 Basic formalism 

The hadronic tensor for deeply inelastic scattering in the case of pure photon exchange is given 
by : 



W^(x,Q 2 ) = — Im T^faq) 



2tt 

i- J d A xe^(P\J,{x)J u (Q)\P) 

2x Ax 2 



-F 2 (x,g 2 ) 



(i) 



Here x = Q 2 / (2p.q) denotes the Bjorken variable, P is the nucleon momentum, q denotes the 
space-like momentum transfer from the leptons to the nucleon with q 2 = —Q 2 and the Compton 



1 Cf. Ref. [3,4] for recent comparisons of the measured value of a s (M\ ). 
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amplitude reads 

T, u (p,q) = i J d A ze^{P\T[J,{z)JM]\P) ■ 
We now consider the light-cone expansion of the (forward) Compton amplitude [17] 

i / d^T[J Ul (z)J U2 (0)]\ m = ^ 



(2) 



jv\vz 



9viHl9 H2V2Q 9viHlQv2Q(J,2 9u2H2QviQ.Hl ~l~ 9viU2Q.fJ.lQH2 

f Q 2 



■;HN} 



(0) 



+ higher twists . (3) 

The twist-2 contributions to the local non-singlet operators O NS (0) are linear combinations out 
of the operators 

a,Hi,-,HN = ^aV^^-'-L^V, a = 1, 2, . . . , (iVf. — 1) , (4) 

/i 2 denotes the factorization scale and a s = a s /(4ir) = g 2 / (An) 2 , with g the strong coupling 
constant. The functions C" L 2 ) N are the Wilson coefficients associated to the moment index N. 
The Mellin moments of the structure functions F 2 N f (x, Q 2 ) are given by 



Mi 



k,N-2 



dx x N 2 F^ s (x, Q 2 ) — C^ S N ( ^ 2 , a s J ^nuci,iV' 



k = 2,L . 



(5) 



For pure electromagnetic interactions is even in x. The corresponding crossing relation has 
the consequence that only the even moments contribute to (3). The nucleon matrix elements 
are 



/ p2\ 

(P|0 NS ' {^ 1 '---'^ JV J'|p} = p{^i---p^iv}^4 NS i N ( j 



(6) 



The scale—dependence of the non-singlet coefficient function is governed by the renormalization 
group : 



r m (91 n 



o, 



2,L . 



The /9-function rules the scale dependence of the coupling constant a s (ji 2 ) 

>da s (^ 



= 0(«O = -£A«i + V) 



(7) 



(8) 



d » 1=0 

on the renormalization scale \x = /xr, where are the expansion coefficients. To 3-loop order 
the constants A for SU(N) [18] 

Po 

Pi : 

#2 = 



34 



3 

2857 



Ci - AC F T F N F - ^C A T F N F 



20 

y 



205 



C\ + 2C 2 T F N F - Z^L Cf C a T f N f 
54 9 



' !15 CiT F iV F + ^C F T 2 N 2 + ^C A T 2 N 2 



27 



158 



(9) 
(10) 

(11) 



contribute, with C A = N C ,C F = (iV c 2 - 1)/(2N C ), and Tp — 1/2 and N c = 3. N F denotes the 
number of flavors. In the QCD-improved parton model the forward Compton amplitude (2), 
valid for nucleon states \P), reduces to that for photon-quark scattering, since only one initial 
state quark participates in the scattering process, with p the quark 4-momentum, 

T<™ = i J d 4 z (p\T [Uz)J v QS)\ \p) ■ ( 12 ) 
One applies the projector 

■ {/il.../ijv} qn 



P 



N 



N\ dp^ ■ ■ ■ dp» N 
onto the Mellin moments iV and the Lorentz projectors 



(13) 

J p 2 =0 



,2 



= (14) 

valid in D = 4 — 2e dimensions. These projections lead to the following moments of the Compton 
amplitude 

£ ) = £ ) ^ 7) k — 2,L. (16) 

3 Renormalization 

The coefficient functions C and Z-factors in (16) obey the following representations : 

C(a , e) = 5 + a (c w + eC u + £ 2 d 2 ) + a 2 , (C20 + eC 2 i) + ajjCso + 0(4) , (17) 

rvl \ 1 , ^11 , 2 /^22 ^21 \ 3 /^33 . ^32 ^31 \ . ~/ 4\ / nQ \ 

Z(a,e) = l + a + a — + + a — + — + + O(a ) , (18) 



V e 2 e J u V e 3 e 2 

with 5 = 1 for C2 and 5 = for Cx,. Here, ao denotes the bare coupling constant, which is related 
to the running coupling by 

a = a - — cf + I -j - — I a + 0(a ) , (19) 

cf. (8). We identified the scales fx 2 = Q 2 . Yet a separation between the contributions to the 
anomalous dimension and the coefficient functions is possible as outlined in the following. The 
anomalous dimension and the Wilson coefficient are 

00 

7 (a,iV) = Y> fc+1 7fc (20) 

k=0 

00 

c(a,N) = 5 + Y,a k c k0 . (21) 
fc=i 
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The respective coefficients 7^ and Cko are determined as follows. We denote by C(£) the coefficient 
of type £ in T (16). Identifying the corresponding powers in e one obtains to 0(a 3 ) : 



with 



70 = 

71 = 

72 = 
ClO = 

C20 = 

C30 = 



Furthermore the relations 



CI? 



c a - 



c a - 



70C10 

C10Z21 — CIIZ22 — 70C20 
C (as ) 

C (a 2 e°) - 7o C n 

C (a 3 e°) — C11Z21 — C12Z22 — 70C21, 



^21 

■^22 

C n 
C12 
C21 



C 



- 7oCio 



C7(ae) 

C (a 2 £ ) - 70C12 



Z3I + ZuC2Q + ^22^11 + Z21C10 



Z32 + Cio^: 



22 








(22) 
(23) 

(24) 

(25) 
(26) 
(27) 

(28) 

(29) 

(30) 
(31) 

(32) 



(33) 
(34) 



hold. The above relations yield the anomalous dimensions and the Wilson coefficients to 0(a 3 ). 



4 Results 

The 16th moment of the 3-loop non-singlet anomalous dimension 7ns + > which describes the 
evolution of the combination 

g+ s (x,g 2 ) = ( gj (x,Q 2 ) + q 3 (x,Q 2 )) - (q k (x,Q 2 ) +q k {x ) Q 2 )) (35) 

of quark densities, and the coefficient functions C^f 16 were calculated using the MINCER algorithm 
[14]. The calculation was performed majorly using two dual-processor 32bit PC's (3 and 2.6 
GHz). Several diagrams were run on an Opteron 64bit PC. Due to the large disk-space required 
by a series of diagrams a 4.2 Tbyte raid system was linked to the two 32bit PC's to store 
intermediary results. In the calculation the moment (16) is determined in terms of an e-expansion 
for a s = a . The anomalous dimension and moments of the coefficient functions are determined 
as described in the previous section. 

Up to 3-loop order 388 diagrams contribute effectively, if symmetry relations between dia- 
grams are used. The calculation using the above system needed about 560 CPU days. The CPU 
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time distribution over the individual diagrams is shown in Figure 1. Several diagrams required 
computation times of a month, in one case of O(60) days. Despite of these long computation 
times the system ran completely stable. The use of a parallel facility running FORM [19] would 
be highly desirable for computations of similar size in the future. 

The 16th moments of the non-singlet anomalous dimension and coefficient functions for 
unpolarized nucleons are : 

, 6(0)i+ = 64419601 = Q207507i 

,NS 6126120 v ; 

i 6)( i )i+ 21546159166129889 3689024452928781382877 ^ 2 

7ns ~ 484994628518400 F A ~ 459818557352009856000 F 

1176525373840303^ Ar 

C F N F 

112588038763200 

= 163.4395247- 13.93310085 N F (37) 

i 6)(2))+ /58552930270652300886778705063429867 59290512768143 \ 3 

7ns ~ V 3451337970612452534317096673280000 ~~ 1563722760600 ^7 F 

( 1670423728083984207878825467 59290512768143 \ 2 
+ V 6488959481351563087872000 + 3127445521200 ^) F A 

/ 1229794646000775781127856064477 59290512768143 \ 2 
+ V 30335885575318557435801600000 ~~ 1042481840400 F A 

( 71543599677985155342551355451 64419601 _ \ ^ _ 

I / I /■ j f^A AT 

V 938967886855098206346240000 765765 S V F F 

( 15018421824060388659436559 64419601 A \ ^ „ Ar 

+ C3 C F C A N F 

\ 579371382263532418560000 765765 s j 

5559466349834573157251 ^ » 

2069183508084044352000 F F 

= 2849.5632736921273714 - 463.86001156801831223 N F 

-3.5823897546153993659 N F . (38) 



£<NS,16 



, , , 4047739719 ^ 

(a s ) = H C F a s 

190590400 

/ 44426674163044428879366970127 24439538 

+ 



+ 



+ 



+ 



V 321931846921747956461568000 ' 255255 
'17918308408498294222783087 113298677 



\ 59422705873182812160000 
143568372761907472111177 



1021020 



C3 )C F Ca 



C F N F 



2758911344112059136000 

/3036813397599509725084677293842505976559161689 
V 8034458016040775933421647863403347968000000 



6 



1494341926940450865387403 59290512768143 27643576 \ 3 
+ 595674040206012768000 ^ 3 + 3127445521200 ^ 4 21879 F 

/262865377883475726558800935515033190333 15355050469171482313 
+ V 56646805852503848671021043712000000 4991403051835200 ^° 

59290512768143 A 47187263 A \ ^ ^ 2 

6254891042400 51051 / ^ ^ 



+ 



+ 



/7750026627118768752845091760890051465242741 
V 1652500620329242273431025887166464000000 
2849482004138921491531 . 59290512768143 . 983963 , \ ^ „ 

6741167121672984000 2084963680800 21879 / F A 

( 4073207241348493196152222079933557529 552298563960959 . 



3529777469944553728278848870400000 4021001384400 



64419601 \ 
+ "1531530" C4 ) ° fNf 

( 582811634921542995647179358698536547 598788865585667 
+ ^ 404620041803598919078721740800000 + 1850495446800 Cs 

1531530 S y 

/ 7227384935999670312318789884999 64419601 



76056398835262954714045440000 20675655 



Ca }C F N 2 F 



/ 705894258514655486993 38404365803 14560 \ d 2 ab 
+ ^ V 3248429831350704000 + 1533061530 ^ 3 5T~ ^7 F 



1 + 28.31719904 a s + (1122.549565 - 69.38406971 N F ) a 2 
+ (50309.36422 - 6651.875513 N F + 131.6959033 N 2 F 

-216.0757466 (e)N F ) a 3 s (39) 



£yNS,16 



+ 



+ 



+ 



+ 



4 

17 



— C F a s 



29393927457809 96 , , 
44659922042736 + Tr 3 ' F + 



/ 55969347000169 
V 8209544493150 



48 

17 Cs 



C F Ca 



39366889 



C F N F 



a. 



39054015 

7508281821276771498126447290110919 
13647898235438852429242598400000 



196256899828170631 
133698296031300 



39360 
17 



/296045501010133565322039207159677 2253147763389895 



V 936620467137960460830374400000 



1188429298056 



Ca 



7 



( 1460792499427100139493280371 1634895686765221 
+ ^ 8256042197255336964480000 2673965920626 ^ 

+ 17 C^jC 2 A C F 

/3529137346321170453160463 _ 44651224 \ 2 
+ V 136796020812222932160000 ~ 765765 ^) F F 

( 4495805144658565385501573689 43594330672 \ „ „ Ar 

+ 1 C3 C A C F N F 

\ 57792295380787358751360000 1249937325 J 

895967716232 ^ - 
209134250325 F F 

, v / 1798450729620489619601 28854977192 2560 , \ d 2 nhr _ 

+ e) (3 H C5 —N f 

w V 18272417801347710000 547521975 S 17 S J N c 

= 0.3137254902 a s + (24.59183098 - 1.344015086 N F ) a 2 s 

+ (1779.007454 - 222.2146165 N F + 5.712233265 N 2 F 

-24.99942737 (e)N F ) a 3 s , (40) 



with 

o N F 

(e) = F E *f (41) 

and d 2 abc /N c = 40/9 for SU(3) C . These results agree with the complete 3-loop results for the 
anomalous dimension in [12]. The moments for the iV^-terms [20] and iVi?-terms [21] were 
known before. For the coefficient functions the moments agree with the very recent complete 
results [15] and an upcoming paper [16]. 



5 Conclusions 

We calculated the hitherto unknown 16th moments for the 3-loop non-singlet anomalous dimen- 
sion 7ns + and the non-singlet coefficient functions C2,l{x,Q 2 ) for pure photon exchange. The 
computation was performed using the MINCER algorithm, which is different from the algorithms 
used in the recent complete calculations. In view of the rather long CPU time of about 560 
days spent for the calculation the reliability of the formula manipulation program FORM has been 
tested intensely as a by-product. The results agree with recent and upcoming complete results. 
The computation of the 16th moment provides a non-trivial test of these computations. 
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Diagram 



Figure 1: Execution time profile for the set of diagrams. Full line: g^ v projection; dash-dotted line: 
P^Pv projection. 
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